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Abstract 

Using data from the E791 fixed-target hadroproduction experiment at Fermilab we have studied 
the Cabibbo-favored but phase-space suppressed decay, D° — ► K~K~K + ir + . We find the decay 
rate for this mode to be (0.54 ± 0.16 ± 0.08) x 10~ 2 times that for the normalization mode D° — > 
K~ir~ir + ir + . We observe a clear signal for D° — ► 0K~7r + which is consistent with producing 0.7 ± 
0.3 of the _D° — > K~K~ K + n + signal. In the context of simple models, we use our measurements 
to estimate the importance of decay amplitudes that produce extra quark-antiquark pairs from the 
vacuum relative to those that do not. 

PACS numbers: 13.25.Ft 
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The decays D° — > K~K~K + tc + and -D — > i^~7r~7r + 7r + are both Cabibbo-favored, but 
phase-space suppresses the former relative to the latter. In addition, the decay D° — > 
K~K~K + n + requires the production of at least one extra quark- ant iquark pair, an ss, either 
from the vacuum or via a final state interaction. The more common decay D° — > K~7t~7t + tt + 
may proceed both via an intermediate state such as K*°p° in which the resonant particles 
contain only quarks produced directly from a spectator amplitude, and via an amplitude 
that requires the production of at least one extra qq pair from the vacuum. 

In this paper we present a decay rate measurement for D° — > K~ K~ K + tt + relative to that 
for D° — > K~7t~7t + 7t + using data from the E791 fixed-target hadropro duct ion experiment 
at Fermilab. This allows us to determine the importance of decay amplitudes that produce 
extra qq pairs from the vacuum relative to those that do not. In addition, we study the 
K~K + invariant mass distribution in signal events to search for intermediate <ft production. 
Experimental Overview 

Experiment E791 recorded 2 x 10 10 interactions during the 1991/92 fixed-target run 
at Fermilab using a 500 GeV/c tt~ beam and an open geometry spectrometer |l| in the 
Tagged Photon Laboratory. The target consisted of one platinum foil and four diamond 
foils, separated by gaps of 1.34 to 1.39 cm. Each foil was approximately 0.4% of a pion 
interaction length thick (0.5 mm for platinum and 1.6 mm for carbon). The average decay 
length of an 80 GeV D° is approximately 5 mm, so most of the D 0, s decayed in the air 
gaps between target foils where backgrounds are lower. Six planes of silicon microstrip 
detectors (SMDs) and eight proportional wire chambers (PWCs) were used to track the beam 
particles. The downstream detector consisted of 17 planes of SMDs for vertex detection, 
35 drift chamber planes, two PWCs, two magnets for momentum analysis (both bending 
in the same direction), two multicell threshold Cerenkov counters[2] for charged particle 
identification (with nominal pion thresholds of 6 GeV/c and 11 GeV/c), electromagnetic and 
hadronic calorimeters for electron identification and for online triggering, and two planes of 
muon scintillators. An interaction trigger required a beam particle and an interaction in 
the target. A very loose transverse energy trigger, based on the energy deposited in the 
calorimeters, and a fast data acquisition system Q allowed the experiment to collect data 
at a rate of 30 Mbytes/s with 50 /is/event dead time and to write data to tape at a rate of 
10 Mbytes/s. 
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Event Selection 

Data reconstruction and additional event selection were done using offline parallel pro- 
cessing systems |4j . Events with evidence of well-separated production (primary) and decay 
(secondary) vertices were retained for further analysis. Candidate D° — > K~K~K + tt + 
and D° — ► fr~7r~7r + 7r + decays (and their corresponding charge conjugate decays, which 
we include implicitly whenever we refer to a decay chain) were selected from events with 
at least one candidate four-prong secondary vertex. Selection criteria (cuts), used for 
both modes, were determined by optimizing the expected statistical significance of the 
D° — > K~K~K + 7i + signal. To avoid bias, we masked the signal region (1.845 GeV/c 2 
< m&ss(K~ K~ K + 7T + ) < 1.885 GeV/c 2 ) in the real data early in the analysis and systemat- 
ically studied sensitivity using a combination of real data for background and Monte Carlo 
simulations for signal. Only after we had determined a set of cuts for the final analysis, and 
looser and tighter sets of cuts for studies of systematic uncertainties, did we examine the 
data in the signal region. 

We used Monte Carlo simulations of D° — ► K~ K~ K + -k + decays and D° — ► K~7t~tt + 7t + 
real data to estimate the efficiencies of potential cuts. We compared Monte Carlo simulations 
of D° —>■ K~7t~7t + 7t + to the real data in this channel to validate our Monte Carlo. Where 
the distributions in the Monte Carlo simulation of D° — > K~tt~7t + 7t + match the real data, 
we trust the D° K~K~K + tt + Monte Carlo simulation. Where the distributions in the 
Monte Carlo simulation of D° — ► K~K~K + tt + match those in the Monte Carlo simulation 
of D° — > K~7t~7t + 7t + , we trust that the corresponding distributions observed for real D° — > 
K~7t~7t + 7t + correctly predict those for D° — > K~K~K + ir + . At this stage, we estimated 
the background within the signal region by interpolating linearly from sidebands in the 
K~K~K + ir + invariant mass distribution of the data. 

When events were initially reconstructed, a topological vertexing algorithm was used to 
identify a primary (interaction) vertex and possible secondary (downstream) vertices. Our 
D° candidates are constructed from four-track secondary vertices (referred to as 'SEED4' 
candidates) and from three-track secondary vertices with the addition of a fourth track 
(referred to as 'SEED3' candidates). Because the initial topological vertexing algorithm 
assigned each track to one vertex candidate, and because it was optimized for two-body and 
three-body charm decays, about half of our signal comes from the SEED4 sample and half 
from the SEED3 sample. 
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FIG. 1: KKKir invariant mass distributions 
after the preliminary event selection 

The signal was expected to be small (between 10 and 50 signal events, depending on 
cuts and allowing for some uncertainty in the relative branching ratio), and the signal-to- 
background ratio is better in the SEED4 sample when all other cuts are fixed, so the cuts 
are studied separately for SEED4 and SEED3 candidates. Within each category, we made 
initial, very loose cuts on a series of candidate variables, primarily informed by our prior 
experience with similar analyses. The invariant mass distribution of all candidates surviving 
these loose cuts is shown in Fig. ^ At this point in the analysis we decided to proceed 
blindly: to mask the signal region in the data and optimize sensitivity as discussed in detail 
below. 

As part of the optimization, we create combined figure-of-merit (FOM) variables from 
kinematic and particle identification variables that are essentially independent of each other 
in our Monte Carlo simulations and in the D° — > K~7i^7T + 7i + data. For each of the SEED4 
and SEED3 samples, we start with a set of independent cuts on individual variables, and 
then construct combined FOM variables for each of the surviving events. The procedure 
for constructing FOM variables will be described in detail after a discussion of the most 
important variables used in selecting candidates. 

The decay vertex is required to lie outside the target foils and other solid material and to 
be well-separated from the primary vertex, with Az > 10 a z for SEED4 candidates (where 
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a z is the error on Az, the longitudinal separation between two vertices) and with Az > 12 a z 
for SEED3 candidates. The transverse momentum of the D° candidate with respect to the 
line-of-flight defined by the secondary and primary vertex positions (p^-balance) is required 
to be less than 250 MeV/c. The transverse distance of closest approach of the D 0, s line- 
of-flight with respect to the primary vertex point (denoted DIP) is required to be less than 
60 jttm. Because vertex separation, transverse momentum imbalance, and DIP correlate 
strongly, we use a very loose DIP cut at this stage of the analysis and include only DIP, of 
these variables, in FOM. 

The sum of the squares of the momenta of the individual tracks with respect to the 
D° momentum vector discriminates between signal and background when normalized to 
the maximum value possible for a candidate's reconstructed mass. (This normalization is 
required to avoid kinematic biases that could artificially create a signal by preferentially 
increasing the acceptance in the signal region relative to acceptance nearby in mass.) This 
ratio is required to be greater than 0.2 for the SEED4 candidates and greater than 0.3 for the 
SEED3 candidates. The product of the ratios of the four daughter tracks' transverse sepa- 
rations from the secondary vertex relative to their transverse separations from the primary 
vertex is required to be less than 0.005 for SEED4 and SEED3 candidates. The maximum 
ratio for a single track is required to be less than 1.0 and is included in FOM for SEED3 
candidates; this is unnecessary for SEED4 candidates because their distribution was cleaner 
upon initial selection. 

To avoid problems due to congestion near the primary vertex, we also found it useful to 
require either an absolute separation of the primary and secondary vertex candidates or that 
the secondary vertex be "isolated" from other tracks by requiring that all other tracks pass at 
least a minimum distance from the secondary vertex candidate. The Monte Carlo simulation 
fails to describe the distribution of additional tracks in the events sufficiently well, so we 
base these cuts, and the associated efficiencies, on studies of the real D° — > K~it~it + it + 
data. For SEED4 candidates we require either vertex separation greater than 0.5 cm or 
secondary vertex isolation greater than 20 /im. For SEED3 candidates we require either 
vertex separation greater than 0.5 cm or secondary vertex isolation greater than 80 /im. 

Our particle identification algorithm compared the light observed in the two multicell 
threshold Cerenkov detectors with that expected for the five hypotheses e, /x, n, K, and 
p for each track. It then assigned a probability to each hypothesis, including a priori 
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likelihoods for each species, so that the sum of probabilities for each track added to unity. 
Tracks that are n/K ambiguous had K probabilities near 0.13. Similarly, those that are 
K/p ambiguous had K probabilities near 0.80. Each kaon candidate is required to have 
K probability greater than 0.20. The three K probabilities are also included in FOM as 
independent variables. Our Monte Carlo simulation of the Cerenkov identification does not 
match our data's variation with transverse momentum well; furthermore, the Monte Carlo 
distributions of K probability for the kaons in D° — > K~K~K + n + and D° — > K~n~n + n + 
decays differ. In calculating the Cerenkov probabilities' contributions to FOM, for each range 
of K probability we use the average of the fraction found in the D° — > K~ K~ K + tt + Monte 
Carlo simulation and in the real data's D° — * K~7t~7t + tt + signal. We also considered using 
the product of the three K probabilities in FOM, but found that the greatest sensitivity 
could be achieved using them independently. Because most tracks are pions, we found no 
benefit in using Cerenkov identification for the pions. 

Two more variables provide some discriminating power between signal and background. 
In the Monte Carlo simulation, there were no SEED4 signal events with proper decay time 
greater than 3.5 ps; in the data, background was observed in this region. So we removed 
SEED4 events with proper decay time greater than 3.5 ps. SEED3 signal events extended 
past 3.5 ps, but we could not find a cut that improved the sensitivity. For both SEED4 and 
SEED3 events, we found that the distribution of the cosine of the polar angle of the sphericity 
axis of the candidate relative to the candidate's line-of-flight discriminated between signal 
and background, primarily because background accumulates preferentially at values of the 
cosine near one. This correlates with the scaled, summed transverse momentum squared 
being small. Having made an absolute cut on the latter quantity, the cosine of the sphericity 
angle is included in FOM for both the SEED4 and SEED3 samples. 

To create FOM variables, we divided the distribution for any one variable into four or 
five ranges and determined the fraction of signal that appeared in each range. Similarly, we 
determined from the data what fraction of background appeared in each range. We use Sa,% 
to denote the probability that a signal event falls in range i for variable A. We use Ba,i for 
background similarly. For example, if three variables A, B, and C are used to define FOM, 
and they are observed in ranges i, j, and k respectively, we calculate: 

FOM = f A '^ BjS c ' k . (1) 

"A,iOB,jOc,k 
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FIG. 2: Figure-of-merit (FOM) distributions for Monte Carlo signal (solid lines) and for 
background (dashed lines) for SEED4 candidates (on the left) and SEED3 candidates (on 
the right), after all cuts on individual variables, as discussed in the text. 

This is effectively a ratio of likelihoods constructed as a product of independent relative 
probabilities. Assuming the variables A, B, and C are statistically independent, an event's 
FOM is the relative signal-to-background ratio in the k cell of A, B, C phase-space: the 
signal-to-background ratio in that cell will be the signal-to-background ratio before FOM 
cuts times FOM. In selecting events using FOM, we accept events with FOM greater than 
a given value and reject those with lower values. 

The FOM distributions for SEED4 and SEED3 Monte Carlo signal and for real data in the 
sidebands (after all the non-FOM cuts) are shown in Fig. El The background accumulates 
preferentially at lower values of FOM, while the signal populates the FOM distribution much 
more uniformly. In determining where to cut on FOM, we calculate the expected KKKw 
signal using (i) our observed K~7t~7t + tt + signal, (ii) the ratio of decay rates previously 
reported by E687p|, and (iii) the relative reconstruction efficiencies determined from our 
Monte Carlo simulations. We calculate the background expected in a 15 MeV/c 2 window 
by extrapolating the K~ K~ K + n + rate from outside our masked-off range (1845 - 1885 
MeV/c 2 ). With no FOM cuts, these assumptions predict 6.1 SEED4 and 7.5 SEED3 signal 



9 





Signal 


Mass 


Width 


real 
data 


KKKir 
Kirinr 


18.4 ± 5.3 
14472 ± 134 


1.8639 ± 0.0015 
1.8658 ± 0.0001 


0.0045 ± 0.0014 
0.0100 ± 0.0001 


Monte 
Carlo 


KKKir 
Kiririr 


595 ± 26 
2156 ± 48 


1.8646 ± 0.0002 
1.8644 ± 0.0002 


0.0041 ± 0.0001 
0.0082 ± 0.0002 



TABLE I: Parameters determined by fitting the final real data and Monte Carlo KKKir and 
Kittttt samples as described in the text. Each Monte Carlo sample was generated with 500,000 
events. The errors quoted are statistical only. 

events with 4.0 and 7.2 background events, respectively. Adding these together predicts 
S/B = 1.2, S/yB = 4.1, and S/y/S + B = 2.7. Our final selection of cuts balances our 
interests in maximizing the S/ y/S + B, maximizing S/ y/B, and maintaining good S/B ratios 
in the SEED4 and SEED3 samples. Our final selection of cuts is FOM > 0.5 for SEED4 
candidates and FOM > 1.0 for SEED3 candidates. With these cuts our algorithm predicts 
11.5 signal events and 3.1 background, giving S/B = 3.7, S/y/B = 6.5, and S/ y/S + B = 
3.0. 

Several points deserve emphasis. The technique for selecting cuts is almost unbiased. The 
data in the D° — > K~ K~ K + tt + signal region have not been examined, so we avoid choices 
which are subconsciously chosen to either enhance the signal level or increase the sensitivity 
of an upper limit should no signal be observed. The potential bias in selecting cuts while 
looking at the background outside the signal region is small; this will be quantified when we 
discuss systematic errors. Using FOM to combine variables that discriminate between signal 
and background allows us to create a relatively robust variable, S/ y/S + B, which correlates 
with our ability to measure the decay rate and varies slowly with changes in FOM cuts. 
We can choose looser and tighter cuts for which S/ B will vary substantially, but S/ sJS + B 
should not. This allows us to examine the data a posteriori to identify potential problems 
with the analysis. 

Branching Ratio Measurement 

The K~ K~ K + tt + invariant mass distribution for events satisfying the final set of cuts 
described above is shown in Fig. El and the i^ _ 7r~7r + 7r + invariant mass distribution used for 
normalization is shown in Fig. |3J The cuts used for the normalization sample correspond 
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FIG. 3: The KKKn invariant mass distribu- 
tion for events satisfying the final selection 
criteria. The SEED4 and SEED3 samples 
have been combined. The fitted signal level 
is 18.4 ± 5.3 events. 
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FIG. 4: The Ktttttt invariant mass distribu- 
tion for events satisfying the selection criteria 
described in the text. These criteria are sim- 
ilar to those used for the KKKn candidates 
to reduce the systematic uncertainties in de- 
termining the relative branching ratio. The 
fitted signal level is 14472 ± 134 events. 



closely to those used for the K~ K~ K + n + sample without FOM cuts. The detailed require- 
ments for finding the vertex outside the target foils or other solid material, and an additional 
requirement that the K~n~n + 'K + daughter tracks not point back to the primary vertex, dif- 
fer slightly because the Q-values of the two decays (summed kinetic energies of the decay 
products) differ substantially. Parameters for the K~ K~ K + n + invariant mass distribution 
are determined using an unbinned maximum likelihood fit in which the signal is described 
as a Gaussian with the mass and width allowed to float and the background is described as 
a quadratic function. Parameters for the K~ K~ K + n + Monte Carlo data as well as for the 
K~n~n + n + real and Monte Carlo data are determined using binned maximum likelihood 
fits in which the signals are described as Gaussian distributions with masses and widths 
allowed to float and the backgrounds are described as linear functions. The quadratic term 
found in fitting the K~ K~ K + n + is small, but we allow this extra degree of freedom to be 
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conservative. The Monte Carlo data has essentially no background, and adding a quadratic 
term to the K~tt~7t + 7t + fit makes a negligibly small difference, so for the Monte Carlo and 
normalization samples we present the results of fits with only linear background terms. The 
results of these fits are summarized in Table HI The K~K~K + tt + signal level is 18.4 ± 5.3 
events. 

To convert this signal level into a ratio of decay rates we need the K~tt~7t + tt + signal 
level (14472 ± 134 events) and the relative efficiency for the two decay modes. The latter 
differs from unity for three reasons. First, the Q- value for the KKKtt decay is smaller than 
that for the Ktttttt decay. This leads to very different track opening angles, and hence to 
very different vertex resolutions. As a result, vertex reconstruction efficiencies and vertex 
separation distributions differ. Second, the background in the KKKtt sample is reduced 
using FOM cuts, a procedure not necessary in the normalization sample. Finally, the KKKtt 
sample has two additional kaons, which reduces the particle identification efficiency We start 
with the relative efficiency determined from the Monte Carlo simulations, 0.275 ± 0.013, 
(where the reported error is the statistical uncertainty from the Monte Carlo samples) and 
make corrections to account for differences between real data and Monte Carlo data observed 
using the D° — * K~tt~tt + tt + signal. These corrections are summarized in Table HT1 Taken 
together, we estimate the efficiency for D° — > K~ K~ K + tt + relative to that for D° — > 
K~tt~tt + tt + to be 15% less than that determined directly from the Monte Carlo simulations. 



systematic variation due to 



raise e K KKir rel by 



correction factor 



kaon Cerenkov efficiency 
SEED4 and SEED3 fractions 
vertex separation 
PT-balance 



+2.5% 
+1.1% 
+0.5% 



0.82 
1.025 
1.011 
1.005 



total correction factor 



-15% 



0.85 



TABLE II: Summary of corrections to the relative efficiency from Monte Carlo for reconstructing 
KKKir and Ktttttt final states. The relative efficiency used for determining the relative branching 
ratio, TKKKir/^Kmrn, will be 0.85 times that found from the Monte Carlo. The total correction 
factor has been calculated by multiplying the individual correction factors. 

The significant sources of systematic uncertainty in the ratio of decay rates are sum- 
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marized in Table IIHI Each of the correction factors described above has a corresponding 
uncertainty determined by studying the D° — > fC~7r~7r + 7r + data that pass cuts, and the 
Monte Carlo samples for both decay modes. The Monte Carlo statistics contribute a 4.7% 
uncertainty. Systematic differences in tracking and vertexing efficiencies between real data 
and Monte Carlo samples have been studied previously p . They contribute an additional 5% 
systematic uncertainty beyond that determined for the correction factors; this is effectively 
an uncertainty on the relative efficiency for the loosest cuts used. We vary the fit used to 
extract the number of D° — > K~ K~ K + ir + signal events allowing both linear and quadratic 
backgrounds, fixed and floating masses, and fixed and floating Gaussian width, and all pos- 
sible combinations. The number of signal events ranges from 17.6 ± 4.7 to 18.8 ± 5.2. Each 
fit describes the data adequately - for each fit the \ 2 P er degree of freedom is less than 
one. We investigated backgrounds due to misidentified charm decays using Monte Carlo 
simulations and found the overall shape to agree well with that found in the fits. We could 
associate a systematic uncertainty of 4.3% with our fitting procedure; this would cover the 
largest excursion of the fit results from the central value reported. Because the background 
near the signal region may be higher than predicted by our fit over the whole range shown, 
we considered several piece-wise linear fits as well. We estimate an additional one event 
systematic uncertainty in the background level which we add in quadrature to give a total 
systematic uncertainty of 7.6% due to the shapes of the signal and background. 

With the corrected relative efficiency described above, and adding the systematic uncer- 
tainties in quadrature, the ratio of decay rates is: 

T( fr = (0.54 ±0.16 ±0.08) x 1Q-* . (2) 

The first error is statistical; the second is systematic. Using the D° — ► K~7t~7t + tt + branching 
ratio reported by the Particle Data Group 7], (7.6 ± 0.4)%, and folding its error into our 
final systematic uncertainty, we obtain 

BTZ(D° -> K'K-K + n + ) = (4.1 ± 1.2 ± 0.6) x 10" 4 . (3) 



The ratio of decay rates reported by E687 is (0.28±0.07±0.01)% j3]. The difference between 
our result and the E687 result is (0.26 ± 0.19)%. The fractional errors for the two results 
are 0.32 (this work) and 0.25 (E687) where the statistical and systematic errors have been 
added in quadrature. 
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differences between Monte Carlo 
and real data (after corrections): 

kaon Cerenkov efficiency 
tracking and vertexing efficiencies 
SEED4 and SEED3 fractions 
vertex seapartion requirement 
PT-balance 
sub total 


10% 

5% 

1.8% 

1.1% 

0.5% 

11.4% 


Statistical fluctuations in Monte 
Carlo data, including difference 
between non-resonant and (pKn 


4.7% 


Signal and background shapes 


7.6% 


FOM predictions 


1.5% 


Total (added in quadrature) 


14.6% 



TABLE III: Summary of systematic errors. The total has been calculated by adding the individual 
contributions in quadrature. 

Search for Resonant Substructure 

Two or more of the final state hadrons in a D° -> K-K-K+-K+ final state might be 
the decay products of an intermediate resonance. The signal observed in this experiment 
is small, and the phase-space so small that it will distort the shapes of broad resonances 
that appear as intermediate states. Hence, we have not attempted a coherent amplitude 
analysis similar to the analysis we did for the decay D° — > i^^i^+Tr - ^" 1 " Q| or similar to 



the incoherent amplitude analysis done by MARK III for the decay D° — > K it tt + tt 



■a. 



Rather, we have looked only at K~K + invariant mass distributions of D° — » K~ K~~ K + n + 
candidates (two pairs per candidate). The signal distribution (for events with 1.855 GeV/c 2 
< m(K~K~K + 7T + ) < 1.875 GeV/c 2 ), seen in Fig. |HJ shows an accumulation of entries near 
the mass. In comparison, the corresponding distribution for background events in the 
ranges 1.700 GeV/c 2 < m(KKKn) < 1.845 GeV/c 2 and 1.885 GeV/c 2 < m(KKKn) < 
2.000 GeV/c 2 , seen in Fig. El has a much broader distribution, with only a hint of any 
structure at low mass. 

To understand the nature of our signal better, we generated two Monte Carlo samples. In 
our non-resonant D° — > K~ K~ K + n + simulation, the generated events populate four-body 
phase-space uniformly. In our D° —>■ <f)K~7r + ] <fi — > K~K + simulation, the generated events 
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FIG. 6: K~K + invariant mass distribu- 

FIG. 5: K~K + invariant mass distribution 

tion for non-resonant D° — > K~ K~ K + w + 

for D° — > (j)Kir Monte Carlo events. There 

Monte Carlo events. There are two entries 

are two entries per D° candidate. 

per D° candidate. 

populate the (0, K, it) three-body phase-space uniformly. Both Monte Carlo samples are 
fully simulated and then reconstructed and analyzed as were the real data. The K~K + 
invariant mass distributions corresponding to those for the real data are shown in Figs. El 
and El The two distributions differ qualitatively. 

Without trying to do a real amplitude analysis, we fit the K~K + invariant mass distri- 
bution of Fig. El as an incoherent sum of the shapes of the two Monte Carlo models and of 
the background region. We use a binned maximum likelihood fit in which the two signal 
fractions float freely and the background fraction floats, but we add a term to the likelihood 
function to account for the difference between the background fraction and that determined 
from the earlier fit of the KKKir used for the branching ratio measurement. This fit (which 
is superposed on the data in Fig. |HJ) finds that 0.7 ± 0.3 of the KKKir signal comes from 
<j)Kir decay. This result indicates that intermediate <fr production is an important mechanism 
in D° -> K~K-K + ir + decay. 

Interpretation of the Relative Branching Ratio 

The relative branching ratio determined in Eq. El is small primarily because the Q- value 
of the KKKtc decay is much less than that of the Knnn decay. The phase space for non- 
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resonant four-body KKKir decay, Vl KKKlT) is only 1.43 x 1CT 2 times that for non-resonant 
four-body Kixixix decay, Q Kwn7T . If one assumes that both decays are purely non- resonant, 
the ratio of invariant matrix elements, TZ, is a constant and can be determined using 




(4) 



(e) (f) 
FIG. 7: Feynman diagrams for amplitudes that can contribute to the decays D° — > K~ K~ K + ir + 
and D° — > i<C~7r~7r + 7r + . When one light quark pair (either dd or uu) in the latter is replaced by 
an ss pair, we get the former. 
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Neither decay is purely non-resonant, but we can use the value of 71, determined using 
this equation, to estimate crudely the importance of amplitudes in which at least one qq pair 
is produced from the vacuum. The decay D° — > K~7t~7t + tt + can proceed via amplitudes in 
which the quarks produced in a spectator decay coalesce to form hadrons which, in turn, 
decay strongly to produce four hadrons in the final state. The decay D° — > K~ K~ K + tt + 
cannot proceed via such amplitudes. It requires either that an extra ss pair be produced 
from the vacuum or that long distance final state interactions of hadrons produced at short 
distances produce such pairs. Feynman diagrams for such amplitudes are shown in Fig. [?l 
along with those of corresponding amplitudes for D° — > K~tt~7t + tt + decays where an extra 
uu or dd pair is produced. 

Accounting for the differences in phase spaces using Qkkk-k /^k-ktt-ki and ignoring any 
possible quantum mechanical intereferences, one can write TZ in terms of the probabilities 
that the final states we are considering are produced by amplitudes in which a qq state is 
produced from the vacuum, P m , P d 2, P ss - for uu, dd, and ss respectively, or in which the 
amplitude has no pair produced from the vacuum, P no pa j r . As a first approximation, one 
can imagine a form of isospin symmetry in which P dd = P S s for each amplitude that can 
lead to four charged hadrons in the final state and the corresponding uu amplitude does not 
lead to four charged hadrons in the final state. In this case we calculate 

U = 0.32 = (5) 

P - + P 

1 ss ~ 1 no pair 

in which case 

= 0.47 . (6) 

-'no pair 

If the amplitudes for producing each flavor of qq pair are the same, and the likelihoods 
for producing four charged hadrons in the final state are the same (allowing for resonant 
three-body decays as well as for non- resonant four-body decay), then one might expect P U u 



P dd = P S s- In this case 



P- 

ss _ = 0.90. (7) 



p 

1 no pair 



If the amplitudes with uu and dd pairs produced from the vacuum somehow interfere destruc- 
tively so that the D° — ► K~ti~'k + 'k + decay rate is equal to that which would be produced 
in the absence of these additional amplitudes, 1Z is a direct measurement of 



P- 

ss - = 0.32. (8) 



p 

1 no pair 
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FIG. 8: K K + invariant mass distribu- 
tions for candidates with 1.855 GeV/c 2 < 
m(K-K-K + ir + ) < 1.875 GeV/c 2 . There 
are two entries per D° candidate. The solid 
line histogram is the real data. The dashed 
line histogram is a toy model in which the 
signal fraction is described as 70% from <f)Kir 
and 30% non-resonant, and the background 
scaled from the data, as discussed in the text. 
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FIG. 9: K~K + invariant mass distribu- 
tion for candidates in the background region. 
There are two entries per K~ K~ K + n + can- 
didate. 



A simple measurement of T(D° — > K K K + tt + ) /T(D° — > K n 7r + 7r + ) cannot tell us which 
picture is closest to the truth, although it seems likely that 0.3 < P S s/P n o pair < 0.9. 



Summary 

Using data from Fermilab experiment E791, we have studied the decay D° — > 
K~K~K + 7i + . To avoid bias, the selection criteria for the K~K~K + tt + candidates were 
determined "blindly" - we masked the signal region in the real data and systematically 
studied sensitivity using a combination of real data for background, and Monte Carlo sim- 
ulations and real data in the normalization decay mode D° — > K~7t~7t + tt + for signal. Only 
after we had determined the final set of cuts did we examine the data in the signal re- 
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gion. We observe a signal of 18.4 ± 5.3 events from which we find the ratio of decay rates 
T(D° -> K'K~K + 7v + )/T(D° -> K-n-n + n + ) to be (0.54 ± 0.16 ± 0.08)%. We also have 
examined the K~K + invariant mass distribution of signal events looking for evidence of 
resonant substructure, i.e., D° — > <pK^Ti + ; <fi — > K~K + . Fitting the distribution using an 
incoherent sum of resonant and non-resonant signal shapes plus a background shape, we 
find that 0.7 ± 0.3 of the signal comes from 4>Kit. Finally, using the ratio of non-resonant 
phase-spaces for the two decays as an approximation for the correctly weighted ratio, we 
find the ratio of matrix elements that lead to the signal and normalization final states to be 
1Z = 0.32 ± 0.10. Producing D° — > K~ K~ K + n + requires producing an extra ss pair from 
the vacuum or in a final state interaction. Relating this probability to TZ is highly model- 
dependent, and our measurement does not suffice to distinguish among models. However, 
it seems likely that 0.3 < P S s/P no pair < 0.9. 
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